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J. S. Lintuvuori,* fl D. Marenduzzo, , K. Stratford,' 7 and M. E. Cates" 



Received Xth XXXXXXXXXX 20XX, Accepted Xth XXXXXXXXX 20XX 
First published on the web Xth XXXXXXXXXX 200X 
DOI: 10.1039/b000000x 

7— I ■ 

We propose a hybrid lattice Boltzmann algorithm to simulate the hydrodynamics of colloidal particles inside a liquid crystalline 
C^l host. To validate our algorithm, we study the static and the microrheology of a colloid in a nematic, with tangential anchoring 
O of the director field at the particle surface, and we confirm theories and experiments showing that the drag force in a nematic is 
D markedly anisotropic. We then apply our method to consider the case of a colloid inside a cholesteric, and with normal anchoring 
at the surface. We show that by tuning the ratio between particle size and cholesteric pitch it is possible to control the defect 
C^) configuration around the particle, and to stabilise novel figure-of-eight or highly twisted loops close to the colloid surface. 
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A suspension of colloidal particles in liquid crystals is an in- 
triguing example of soft material, which has been receiving in- 
creasing attention in the past few years among physicists and 
chemists 1 " 5 . Liquid crystals are anisotropic materials whose 
average coarse grained orientation may be described by a vec- 
tor field, known as the director field. The presence of colloidal 
particles creates a disturbance in the director alignment, which 
results in the formation of disclination lines or other topolog- 
ical defects. The interactions between two or more colloids 
embedded in a liquid crystal is governed by the elasticity of the 
surrounding medium and by defect-defect interactions, both 
of which are long range. Experiments, theory and simulations 
have shown that these elastic and topological long range in- 
teractions may provide in practice a route to the formation of 
a striking variety of self-assembled spatially organised struc- 
tures, making colloids in liquid crystals a promising example 
of a material with controlled and programmable macroscopic 
properties. 

For instance, Poulin et al. have shown that large colloidal 
particles in nematics with homeotropic (i.e. normal) anchor- 
ing of the director field at their surface form lines or some- 
times branched networks. More recently, Araki et al. 2 have 
performed simulations of smaller particles, again with nor- 
mal anchoring and in a nematic host, and found that the par- 
ticles form very long-lived aggregates, which are sometimes 
quasi-planar and are held together by a single disclination loop 
embracing all the particles in a cluster. Several theoretical 
works based on liquid crystal elasticity 3-5 have successfully 
explained why particles of different sizes behave so differ- 
ently: this is due to the different defect configuration minimis- 
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ing the liquid crystalline free energy close to the colloid sur- 
face. For small particle size, there is a disclination at the equa- 
torial plane of the colloid, known as a "Saturn ring", which has 
quadrupolar symmetry. The disclination line has topological 
charge -1/2, whereas the total charge of the loop is -1. For 
large size, instead, the minimal energy configuration corre- 
sponds to a point defect of integer topological charge (equal 
to -1) on one side of the particle, so that the resulting colloid 
plus defect system now has a dipolar symmetry. The transi- 
tion between Saturn ring and dipole may also be triggered by 
increasing the strength of the anchoring at a fixed particle size. 

Previous simulation work on the topic has focussed on this 
case of spherical particles in nematics with normal anchor- 
ing 2 ' 5 " 9 . With the exception of the previously referred simu- 
lations in 2 which solve for the hydrodynamics of very viscous 
and rigid droplets, the majority of other numerical work have 
considered a purely elastic model in which the nematodynam- 
ics of the liquid crystalline host is neglected. This method pro- 
vides a way to energy minimisation and as such as has been 
very useful to characterise dimerisation and oligomerisation in 
nematic colloids, as well as the formation of entangled discli- 
nation networks in regular two or three-dimensional arrays of 
particles. On the other hand, the non-hydrodynamic frame- 
work by necessity rules out e.g. rheological studies, which 
are experimentally useful ways to characterise soft materi- 
als macroscopically. Furthermore, it is possible that the self- 
assembly of a large number of colloids may be dynamically 
controlled, in which case the prediction of the final structure 
by theory requires a faithful description of the host hydrody- 
namics. 

Our programme in this work is therefore to take a different 
route to study the hydrodynamics of colloidal particle in liquid 
crystals, by means of hybrid lattice Boltzmann (LB) simula- 
tions. LB has long proved to be a powerful method to study 
colloidal suspensions in Newtonian or binary fluids 10 ~ 14 5 and 



we show here that it promises to be equally effective to probe 
the hydrodynamics of colloids in liquid crystals. Our algo- 
rithm is the first hybrid LB based algorithm to study the hy- 
drodynamics of colloidal particles in liquid crystals. The use 
of an LB formalism to solve the Navier-Stokes equation al- 
lows us to take flow effects into considerations, unlike most of 
the previous work in the literature (with the exception of the 
different algorithm used in 2 which however approximates the 
solid particles with viscous fluid droplets). 

To validate our LB method, we first reproduce earlier ex- 
perimental and theoretical observations that the drag force on 
a sphere moving in a nematic liquid crystal is anisotropic. 
In the second part of our work, we instead apply our algo- 
rithm to find new results on the problem of a colloidal parti- 
cle * with homeotropic anchoring immersed in a cholesteric 
liquid crystal. In particular, we describe the disclination pat- 
tern and the director field surrounding a colloidal particle with 
homeotropic anchoring in a cholesteric host. We show that by 
varying the ratio between particle size and cholesteric pitch 
it is possible to control the disclinations around the particle 
- upon increasing size this changes from a Saturn ring to a 
figure-of-eight to a twisted loop which covers the surface more 
uniformly (note that with the parameters we chose the range of 
particle sizes is such that in a nematic the Saturn ring is always 
the stable configuration for particle with homeotropic (nor- 
mal) surface anchoring). Previous related but distinct work 
had either focussed on colloidal particles in twisted nemat- 
ics 9 , realised via conflicting anchoring at the boundary planes 
and a nematic liquid crystals, or in blue phases 15 , where the 
emphasis was on exploring the interactions between nanopar- 
ticles and the disclination lines making up the blue phase lat- 
tice. Our results instead are relevant to the standard cholesteric 
phase, and our predictions could be tested with microscopy 
experiments in the bulk. 

The study of colloidal intrusions in a cholesteric host re- 
quires the accurate resolution of an additional length scale, the 
cholesteric pitch. In practice, LB simulations of cholesterics 
need several (in our case 16 or more) lattice sites to resolve 
the cholesteric half pitch. As we are interested in the case 
in which the particle size can compete and be larger than the 
cholesteric pitch (this is likely to be the experimentally rele- 
vant case), the colloid needs to be rather large in size. This im- 
mediately leads to numerical limitations as the Brownian time 
scales with the cube of particle size, hence problems in which 
we need colloids to diffuse significant distances are compu- 
tationally costly. These can only be performed via parallel 
simulations in supercomputers. While we are exploring these 
studies, here we limit ourselves to one-particle simulations, 
which do not require such long equilibration times yet are in- 
teresting as they lead to new physics and can already answer a 



* Note that we here use the term "colloid" also to refer to a single particle. 



number of interesting questions. 

2 Simulation method 

2.1 Liquid crystalline equations of motion 

The thermodynamics of the cholesteric solvent is determined 
by the Landau - de Gennes free energy J , or equivalently its 
density /, which is expressed in terms of a (traceless and sym- 
metric) tensorial order parameter 16 Q. Here we briefly review 
the theory, and refer the reader to e.g. 16-18 for more details. 
The free energy density we use is 

/ = £ (1 - \) Gap - ^fi«pfipr0i« + A -f(Qlp) 2 

+ § (Vpgap) 2 + § (ea Y 5V Y e § p +2q Qav) 2 - W 

In the formula above, Aq is a constant, K is an elastic constant, 
go = 2n/p, p being the cholesteric pitch, and y represents 
a temperature-like control parameter for thermo tropic liquid 
crystals. In our notation Greek indices denote Cartesian com- 
ponents and summation over repeated indices is implied. Note 
the limit qo =0 corresponds to a nematic host. The distortion 
free energy in that limit is equivalent to the more commonly 
used j (V a gpy) 2 up to a total divergence. 

In the Beris-Edwards model, the evolution of the Q tensor 
is described by the following equation 16 

AQ = r(^ + lTr(g)l). (2) 

Here T is a collective rotational diffusion constant, |jj is a 
functional derivative and D t is the material derivative for rod- 
like molecules 16 , while Tr denotes the tensorial trace. The 
term in brackets is the molecular field, H, which ensures that 
in the absence of flow Q evolves towards a minimum of the 
free energy. The velocity field, u, obeys the continuity equa- 
tion and the Navier-Stokes equation, 

p (3, + MpVp) u a = VpP a p + r| Vp (V aM p + Vp Ma ) , (3) 

where p is the fluid density and r| is an isotropic viscosity. 
The stress tensor appropriate for liquid crystal hydrodynamics 
is explicitly given by: 

^ap= ft5ap + 2^(e a p + ^S a p)e >E ff ye (4) 

- tflayiQyfi + 5§ 7 p)-t;(ga7+5 §ay)flyp 

— V a Qyv g — h QayHyR — H a yQya , 

° v p£V 

where £, is related to the aspect ratio of the molecules and Pq 
is an isotropic pressure 16 . 

There are a few dimensionless parameters which are useful 
to keep in mind and which determine the physics of the liquid 



crystal host and of the colloidal particle. On one hand, the 
equilibrium thermodynamics of chiral liquid crystals is mainly 
controlled by two quantities, the chirality K and the reduced 
temperature T, which are given in terms of previously defined 
quantities as: 



lOSKql/Aoy 
27(l-y/3)/y. 



(5) 



The physics of a colloidal particle moving in the liquid crys- 
tal is on the other hand primarily controlled by the Ericksen 
number, 



YivR . , 2q 2 
Er = , with Yi = 
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(6) 



which measures the ratio between viscous and elastic forces. 
In the expression for Er, yi is the rotational viscosity of the 
liquid crystal, q is the degree of ordering in the system (for 
uniaxial case with director n, Q a p = q(n a n^, — 8 a p/3)) and 
v and R are the velocity and radius of the colloidal particle 
respectively. 

The Navier-Stokes and order parameter Beris-Edwards 
equations characterising the liquid crystalline host can be 
solved via a 3D hybrid lattice Boltzmann/finite difference 
(LB/FD) algorithm detailed elsewhere 19 . 

In the following simulations we have modelled a liq- 
uid crystal with these parameters in simulation units: Aq = 
1.0, K~ 0.065, £ = 0.7, y= 3.0,4 = 1/2 and T = 0.5. They 
lead to effective temperature T = 0, and rotational viscosity 
Yi = 1. Assuming an elastic constant of 6.5pN ', rotational 
viscosity of 1 Poise (and a colloidal diameter of 1/im), we can 
map the simulations units for force, time and velocity onto 
~ 100 pN, 1/js and 0.03/jm/s, respectively. We note that these 
parameters are close to those of liquid crystals typically used 
in experiments or devices, such as 5CB. 

2.2 BBL and hybrid dynamics 

To represent colloids as moving solid objects within the LB 
method, we have made use of the well-established procedure 
of bounce-back on links (BBL) 10 " 12 . In the hybrid LB/FD 
approach used here, BBL is retained for the LB density dis- 
tributions, which are simply reflected at the solid-fluid surface 
with a correction which depends on the local surface velocity. 
The resulting change in momentum is summed over the links 
to give the net hydrodynamic force on the colloid, which is 
then used to update the particle velocity, and hence position, 
in a molecular dynamics-like step. Boundary conditions for 
the finite difference equations for the order parameter tensor 
are dealt with in a different way. 



t Within the one elastic constant theory used here, the Frank elastic constants 
are given in terms of K as 2Kq 2 . 



First, we note that the assignment of solid and fluid lattice 
nodes for the order parameter follows that for the density: in- 
side and outside are distinguished using the nominal radius 
of the colloid ciq and its position. It is useful, in addition, to 
think about a series of control volumes surrounding each lat- 
tice node whose faces are aligned with the lattice (Fig. 1). 
A set of these faces constitute the solid-fluid boundary in the 
hybrid picture. 

Boundary conditions for Q a p are of two types: 
homeotropic, where the director is aligned with the 
local unit normal to the surface s a , and planar, where the di- 
rector lies in the plane of the tangent to the surface. The exact 
orientation in the planar case is determined by a projection 
of the existing director onto the local tangent plane, with unit 
normal s. (In this case the director is allowed to relax in the 
finite difference step). For either choice of director at the 
surface, we may set the corresponding value of Q a p at lattice 
nodes immediately inside the surface via 



(7) 



where the constant S° controls the degree of surface order. 
This allows us to compute, at all fluid nodes, the derivatives 
Vygap and V"ga|3 using the same finite difference stencil. 
This allows the molecular field and hence the diffusive terms 
in the Beris Edwards equation to be computed. 

Also appearing in the Beris-Edwards equations is the ve- 
locity gradient tensor (this comes in the material derivative 
D t in Eq. 2), which can be handled in a similar fashion close 
to the colloid. The velocity field at solid nodes immediately 
inside the colloid surface are set to the solid body velocity 
u = v + COxr(co denotes the angular velocity of the colloid). 
Again, the velocity gradient tensor 3 a t'p may be computed us- 
ing the same stencil at all fluid nodes. 

Advective fluxes of order parameter are computed at the 
faces of the control volumes, and the boundary condition is 
zero normal flux at solid-fluid interfaces. At the given time 
step, the discrete colloid is effectively stationary, so the no 
normal flux condition assumes a stationary boundary. 

The force on the fluid originating from the order parameter 
is computed via the discrete divergence of the stress P a p. In 
the fluid, this is implemented by interpolating P a p to the con- 
trol volume faces and taking differences between faces in each 
direction. This method has the advantage that, with the intro- 
duction of colloids, an interpolation/extrapolation of P a p to 
the solid-fluid boundary is possible. This allows one to com- 
pute the divergence of the stress at fluid nodes adjacent to the 
colloid. At the same time, the discrete equivalent of 



rcoll _ 
r a ~ 



dS 



(8) 



is found by summing P a p over the relevant solid-fluid control 
volume faces. By construction, this ensures that momentum 



lost by the fluid is gained by the colloid, i.e., global momen- 
tum is conserved. 

Finally, movement of the colloid across the lattice is accom- 
panied by changes in its discrete shape. The events necessitate 
the removal or replacement of fluid information. For the re- 
placement of fluid at newly exposed lattice nodes, this means 
an interpolation of nearby order parameter values in the fluid 
to provide the new information. This is analogous to what is 
done for the LB distributions. 
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Fig. 1 The colloid (represented by the solid circle) moves 
continuously across the lattice. Lattice sites inside are designated 
solid, and those outside fluid (open and closed points, respectively). 
In the lattice Boltzmann picture (left) the surface is defined by a set 
of links fa, which involve discrete vectors Ci,At which connect fluid 
and solid sites. For the order parameter (right), the colloid is 
represented by the set of faces, e.g., that between sites i,j and 
i + 1, j with unit normal S x . Discretisation effects are found to be 
negligible for radii greater than about 5 lattice units. 



3 Results 

3.1 A colloidal particle in a nematic liquid crystal: pla- 
nar anchoring 

To first validate our algorithm, we start from the case of a col- 
loidal particle in a nematic liquid crystal. We consider planar 
anchoring, which has been well studied theoretically and ex- 
perimentally before, but, due to the degeneracy in the allowed 
director configurations, is a more challenging case for sim- 
ulations with respect to the normal anchoring. Indeed most 
simulation work has focussed on the latter 7 " 9 ' 20 (but see Ref. 6 
for an interesting exception using mixed planar and normal 
anchoring on a single "Janus" particle). 

A colloidal particle with planar surface alignment inside a 
nematic liquid crystal leads to a defect structure confined to 
the particle surface, with two surface point defects on oppo- 
site sides of the particle along the far field nematic director 
(Fig. 2(a)*) Each of the defects has integer topological charge 
(equal to +1, in agreement with the Poincare theorem, see 




Fig. 2 (a) A colloid with planar surface anchoring in nematic host 
results in a defect structure with two surface defects of topological 
charge +1 (boojums). (b) 3D plot of the director field close to the 
particle surface with color coding indicating the alignment with the 
far field nematic director, n (marked by the arrow). Blue and red 
denote parallel and perpendicular alignment respectively. Frame (b) 
was produced with QMGA. 21 

e.g. 22 ). Fig. 2(b) shows the director field (individual directors 
are represented by spherocylinders) close to the partice sur- 
face. Here the color coding quantifies the alignment between 
the local and the far field nematic director (marked by the ar- 
row): blue indicates parallel alignment and red perpendicular 
alignment. 

The bounce-back on links allow us to have approximately 
no slip on the particle surface. This is very well achieved ev- 
erywhere apart from close to the boojums, where the max- 
imum spurious velocity is about 0.001 (in LB units) after 
50,000 steps. 

One of the main strength of our LB algorithm is that, on 
top of recovering the correct statical properties of the colloid 
plus director field, it readily allows hydrodynamic studies. We 
exploit this by directly simulating an active microrheological 
experiment, in which a colloid is dragged inside a nematic, 
both perpendicularly to and along the far field nematic direc- 
tor field (again the director is forced to lie tangentially to the 
particle surface). This is interesting in view of the passive 
microrheology experiments and drag calculations for this sys- 
tem, which both predict markedly anisotropic behaviours and 
effective viscosities in the two cases 23 ' 24 . Fig. 3a shows a 
typical time series of the instantaneous velocity of a colloid 
pulled with constant force (in this case / = 0.09) along, and 
perpendicular to, the far field director. Fig. 3b shows the cor- 
responding velocity-force curves. These force-velocity data 
were fitted by the means of the following formula for the drag 
force, 

f = av. (9) 

Here a has the dimension of viscosity, which is different 
for perpendicular and parallel dragging. In the former case, 



f Note that the point defects appear as patches in our rendering as we plot iso- 
surfaces of the order parameter. The defects are then the small volumes en- 



closed by these isosurfaces, within which the order parameter is below a given 
threshold, which in our case was fixed to 0.435. 
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Fig. 3 (a) Plot of the velocity versus time for a colloidal particle 
pulled along (solid line) and perpendicular to (dashed line) the 
nematic director with constant force / = 0.09 (a) (left axis). In the 
same panel the ratio between the two velocities is shown (dot-dashed 
line, right axis). The dotted line is the ratio ax/ a || = T|_l/t| || , 
obtained from the fits in (b) (see the text for details), (b) Plot of the 
force- velocity curve for the system pulled perpendicular (open 
symbols) and along (closed symbols) the nematic director. The fits 
are performed via Eq. (9) in the text. (Particle radius is R = 16, and 
the isotropic viscosity T|, which enters the Navier-Stokes equation, is 
equal to 0.1 in simulation units.) 



a± « 116.9 ±0.2, whereas in the latter an ~ 63.2 ±0.2. As- 
suming an effective Stokes' law to hold for the viscous drag, 
such that / = 6nr[R, we can now interpret our results as point- 
ing to the existence of two distinct effective microrheological 
viscosities r| eff , which can be probed when dragging a spher- 
ical particle inside a nematic liquid crystal. We find that the 
drag force is anisotropic with the ratio between these two ef- 
fective microviscosities perpendicular/parallel to the nematic 
director ri^/ri^ = a , /an w 1.850 ±0.004. This result com- 



: a±/a\\ 

pare quite well to /r|jj ff « 1.6 and /r^ « 1.72 ob- 
served earlier by experiments 24 and predicted theoretically 23 
for a Stokes drag experienced by a Saturn-ring droplet in ne- 
matic liquid crystal at low Er (for simulations with larger Er 
see 25 ' 26 ). However our results hold for planar anchoring in 
which the defect structure is a pair of boojums (Fig. 2) rather 
than a Saturn ring. This points to a rather weak dependence of 
the viscosity ratio on the details of the anchoring of the direc- 
tor fields at the particle surface. 



3.2 A colloidal particle in a cholesteric: normal anchor- 
ing 

As mentioned in the Introduction, by varying particle size 
and/or anchoring strength, one can trigger a switch in the con- 
figuration of the defect close to a particle with normal anchor- 
ing immersed in a nematic liquid crystal, from a Saturn ring to 
a dipolar defect. As discussed several times in the literature on 
colloids in liquid crystals 1 - 9 ' 26 ' 27 ? the defect structure will typ- 
ically determine the interparticle interactions, and hence the 
ultimate self-assembled structure, in a system at finite colloid 
density. Hence it is interesting to ask whether there may be 
other ways to control defect topology by altering other prop- 
erties in a colloidal suspension in a liquid crystal. To explore 
this possibility, we have studied the statics of spherical particle 
of various sizes inside, this time, a cholesteric liquid crystal. 
This case is readily treatable within our framework, by just 
selecting the Landau-de Gennes distortion free energy appro- 
priate for cholesterics. 

Fig. 4 shows the minimal energy configurations found by 
our algorithm. We focus on a range of particle sizes for which 
the free energy inside a nematic host is minimised through a 
Saturn ring arrangement (Fig. 4a). We have then considered 
a cholesteric liquid crystal, with pitch, p = 32, lattice units 
(giving chirality K = 0.3), and various particle sizes. For a 
particle with R = p/4 (so that its diameter is half the pitch), 
the Saturn ring elongates and twists up into a figure-of-eight 
structure. This is similar to a metastable structure found in a 
colloidal dimer in nematics 2 : in this case however this state 
is stable and just requires a single colloid. Increasing the 
size of the particle, the disclination loop twists and becomes 
longer, providing a more uniform coverage of the particle sur- 
face (Figs. 4c,d). 

In order to better compare with existing and future exper- 
iments, it is also useful to plot the cholesteric layers close to 
the particle, as these may be observed with optical microscopy 
experiments. These are shown in Fig. 5, where we have color- 
coded one of the components of the director field in the plane 
perpendicular to the direction of the cholesteric helix (see cap- 
tion of Fig. 5). It can be seen that there is a significant bend- 
ing of the layers right above and below the colloidal parti- 
cle, whereas laterally there is a shrinking and dilation of the 
cholesteric layers. In particular, the layers in which the di- 
rector ordering conflicts with that imposed by the anchoring 
shrink, and this is more marked close to the disclination line. 
What is perhaps surprising is the limited distortion of the liq- 
uid crystal, which extends to less than its size. 

Interestingly, in a recent experimental work performed at 
the Kent State University 28 ^, a colloidal particle was analysed 
while it was sedimenting. It was found that the particle ro- 

j See also the information available online at the Kent State web-site at 
http://www.lci.kent.edU/Lavrentovich/Overview.html#rotation. 




Fig. 4 Defect structure close to the particle surface. In the nematic 
limit (R/p = 0) a Saturn ring defect is formed (a). For cholesteric 
case (R/p > 0) twisted Saturn rings are formed: (b) R/p = 1 /4, (c) 
R/p = 1/2 and (d) R/p = 3/4. (The arrow denotes the direction of 
the cholesteric helix.) 




Fig. 5 Cut of the director field in the yz plane, z being the direction 
of the cholesteric helix, see coordinate axes in the picture. The 
director field is color-coded according to its projection on x (see 
color bar). The disclination lines are also shown in blue. The radius 
of the colloidal particle with is equal to p/2. 



ology of these particles, both in isolation and in concentrated 
solutions. Our LB algorithm is in principle capable of ad- 
dressing these issues, and we hope to explore some of them 
elsewhere. 



tated while translating. More relevant to our analysis, the de- 
fect structure formed around the particle was determined by 
optical microscopy, and it was found that it consisted in spi- 
rals going around the particles, akin to the ones we found with 
our LB algorithm. Furthermore, via fluorescence confocal po- 
larizing microscopy the deformation of the cholesteric layers 
was determined, and the results are in good qualitative agree- 
ment with our snapshot in Fig. 5. It should also be noted that 
our picture refers to a steady state result, for a particle which 
is at rest in equilibrium - it can be imagined that a particle in 
motion such as a sedimenting one may result in larger defor- 
mations. 

We close this section with an overview of possible future 
calculations which our results may open up. In order to un- 
derstand the self-assembly potential of colloidal particles in 
liquid crystals, it would be important to probe the different 
effective interactions that are caused by the defect configura- 
tions which we have described. As mentioned in the Introduc- 
tion, in a nematic liquid crystal Saturn ring forming colloids 
aggregate in clusters kept together by disclination loops. It 
might well be that the twisted chiral Saturn rings found here in 
cholesterics aggregate or interact quite differently. Finally, it 
would also be interesting to simulate the micro and macrorhe- 



4 Conclusions 

In conclusion, we have presented a new hybrid lattice Boltz- 
mann algorithm to simulate the hydrodynamics of colloidal 
particles in a liquid crystal. 

To validate our algorithm, we have first considered the case 
of a colloid in a nematic, with tangential anchoring at the 
director field at the particle surface. We have reproduced 
the known defect structure, consisting of two surface defects 
(boojums) of integer topological charge, and also addressed 
the hydrodynamic problem of the viscous force which is felt 
by one such particle as is dragged inside the liquid crystal. 
As expected from theoretical calculations and experiments fo- 
cussing on Saturn rings, we have found that the drag force is 
anisotropic, with the colloid moving more easily along the far 
field director field than perpendicular to it. The direct counter- 
part of our simulations would be microrheology experiments, 
e.g. performed with optical tweezers in the fixed-force mode, 
where the force felt by the particle is maintained constant via 
a feedback mechanism. 

We have then considered the static of a spherical particle 
embedded with normal anchoring in a cholesteric liquid crys- 
tal, and found that by varying the particle size to cholesteric 
pitch ratio it is possible to control the defect configuration 



close to the surface of the particle. In particular, the Saturn presentation COL25 at the 22nd International Liquid Crystals Confer- 

ring which embraces the particle in the nematic limit twists ence - Je J u - Korea, 2008. 

up and elongates, to form e.g. a figure-of-eight loop when the 
particle size equals the cholesteric half pitch. Furher increase 
in particle size leads to longer and more twisted loops, which 
more uniformly cover the surface. Therefore the defect con- 
figuration may be controlled by tuning the size to pitch ratio, 
and it would be interesting to understand what consequences 
this tunable topology has on the self-assembly of a solution of 
such colloidal particles. 

We acknowledge funding by EPSRC Grants EP/E030173 
and EP/E045316. MEC is funded by the Royal Society. 



References 

1 P. Poulin, H. Stark, T. C. Lubensky and D. A. Weitz, Science, 1997, 275, 
1770. 

2 A. Araki and H. Tanaka, Phys. Rev. Lett, 2006, 97, 127801. 

3 H. Stark, Eur. Phys. J. B, 1999, 10, 31 1. 

4 R. W. Ruhwandl and E. M. Terentjev, Phys. Rev. E, 1997, 56, 5561. 

5 J. Fukuda, M. Yoneya and H. Yokoyama, Eur. Phys. J. E, 2004, 13, 87. 

6 M. Conradi, M. Ravnik, M. Bele, M. Zorko, S. Zumer and I. Musevic, 
Soft Matter, 2009, 5, 3905. 

7 M. Ravnik, M. Skarabot, S. Zumer, U. Tkalec, I. Poberaj, D. Babic, 
N. Osterman and I. Musevic, Phys. Rev. Lett, 2007, 99, 247801. 

8 M. Skarabot, M. Ravnik, S. Zumer, U. Tkalec, I. Poberaj, D. Babic and 
I. Musevic, Phys. Rev. E, 2008, 77, 061706. 

9 U. Tkalec, M. Ravnik, S. Zumer and I. Musevic, Phys. Rev. Lett, 2009, 
103, 127801. 

10 A. J. C. Ladd, J. FluidMech., 1994, 271, 285. 

11 A. J. C. Ladd, J. Fluid Mech. , 1994, 271, 311. 

12 N.-Q. Nguyen and A. J. C. Ladd, Phys. Rev. E, 2002, 66, 046708. 

13 M. E. Cates, K. Stratford, R. Adhikari, P. Stansell, J. C. Desplat, I. Pago- 
nabarraga and A. J. Wagner, J. Phys. Cond. Matt, 2004, 16, 3903. 

14 K. Stratford, R. Adhikari, I. Pagonabarraga and M. E. Cates, Science, 
2005, 309, 2198. 

15 M. Ravnik, G. P. Alexander, J. M. Yeomans and S. Zumer, Faraday Dis- 
cussions, 2010, 144, 159. 

16 A. N. Beris and B. J. Edwards, Thermodynamics of Flowing Systems, 
Oxford University Press, Oxford, 1994. 

17 O. Henrich, D. Marenduzzo, K. Stratford and M. E. Cates, Comput. Math. 
Appl, 2010, 59, 2360. 

18 O. Henrich, K. Stratford, D. Marenduzzo and M. E. Cates, Proc. Natl. 
Acad. Sci. USA, 2010, 107, 13212. 

19 M. E. Cates, O. Henrich, D. Marenduzzo and K. Stratford, Soft Matter, 
2009, 5, 3791. 

20 M. Ravnik and S. Zumer, Soft Matter, 2009, 5, 269. 

21 A. T. Gabriel, T. Meyer and G. Germano, /. Chem. Theory Comput., 2008, 
4, 468. 

22 M. Kleman and O. Lavrentovich, Phil. Mag., 2006, 86, 41 17. 

23 H. Stark and D. Ventzki, Phys. Rev. E, 200 1 , 64, 03 1 7 1 1 . 

24 J. C. Loudet, P. Hanusse and P. Poulin, Science, 2004, 306, 1525. 

25 A. Araki and H. Tanaka, /. Phys. Condens. Matt., 2006, 18, L193. 

26 B. T. Gettelfinger, J. A. Moreno-Razo, G. M. Koenig, J. P. Hernandez- 
Ortiz, N. L. Abbott and J. J. de Pablo, Soft Matter, 2010, 6, 896. 

27 I. I. Smalyukh, O. D. Lavrentovich, A. N. Kuzmin, A. V. Kachynski and 
P. N. Prasad, Phys. Rev Lett, 2005, 95, 157801. 

28 B. I. Senyuk, O. D. Lavrentovich, C. J. Smith and C. Denniston, Poster 



